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ABSTRACT 


The mechanics of separation of a thin interfacial liquid layer trapped between two 
parallel surfaces was studied in a controlled manner. Different liquid viscosities, layer 
thicknesses and separation velocities were used to investigate the extensional behavior 
and determine its dependence on viscous fingering, and Capillary number. Force, 
displacement and time data have been recorded for all experimental runs. Qualitative 
visual data have also been recorded to corroborate the trends in the onset of viscous 
fingering based on a simple interfacial stability analysis. The quantitative data has been 
used to generate force-displacement plots of the separation. The results of this work 


provide useful fundamental insight into the mechanics of this novel problem. 
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I. INTRODUCTION 


The original studies on the instability of the interface between two liquids were 
carried out in the late 1950s, when water was pumped down into an apparently ‘dry’ oil 
field to push up the remaining oil. The belief was that the oil being lighter than water 
would rise, although in reality the interface between the water and oil was unstable and 
fingers of water penetrated up through the ground leaving the oil down below. It was 
shown by Taylor (1958), and Saffman and Taylor (1958), in a now classic study that the 
phenomenon could be explained by a simple analysis based on Darcy’s law for flow of a 
viscous fluid through a homogeneous porous medium. 

The study of the motion of fingers of a less viscous fluid through a more viscous 
fluid in a Hele-Shaw cell (i.e. in a narrow gap between parallel plates.) is a paradigm in 
applied mathematics. It has raised questions at all levels. By choosing the proper 
equations of motion, it can be solved analytically or numerically. Considerable 
experimental evidence has also been gathered for various scenarios. Continuing 
disagreement with experiments leads to modifications of the assumptions and 
approximations employed by the theory until good agreement is reached to help 
understand the phenomenon (Saffman, 1991). 

It has been shown that for a channel of given aspect ratio (width/thickness of the 
layer), it was found for the case when the viscosity of the displacing fluid is negligible 
that the shape and width of the fingers depend only on the Capillary number (Ca). A way 
to measure the ratio of viscous forces to surface tension is through the use of the 
Capillary number (Ca = pU/o). If the capillary number remains small (107°-10°), the 


surface tension effects are expected to dominate the flow. The combination of the liquid 


viscosity and interface velocity could be chosen in such a way that the capillary number 
would be in the above limits. Also, by keeping the aspect ratio large (D/lp >>1), the fluid 
mechanics in the layer could be restricted to essentially one-dimensional behavior. 

In this experiment, the approach to the problem is different. A tensile force is 
applied to a layer of liquid trapped between two circular parallel disks and the force- 
displacement behavior is studied vis-a-vis the onset of viscous fingering. A wide range 
of viscosities, separation velocities and initial liquid layer thicknesses are explored. 

A variation of the Saffman-Taylor instability analysis will be used to determine 
the neutral stability curve and the most dangerous wavelength that grows into fingers, as 
well as the number of corrugations (fingers) that form. Through experimental analysis of 
the data, the interface velocity and the Capillary number at the point of viscous fingering 
can be calculated with extreme accuracy and repetition. 

Quantitative force-displacement data and qualitative visual data have been 
gathered for this experimental study. Both forms of data have been used to make some 


useful deductions of the mechanics of separation of thin interfacial liquid layers. 


Ul. BACKGROUND 


A. SURFACE TENSION 

hk Theoretical Background 

The effects of surface and interfacial tensions give rise to so many phenomena 
observed in liquid behavior that the complex physical-chemical interactions involved are 
taken for granted and not all of which are understood even today. The liquid state itself is 
composed of molecules in motion that are kept relatively close to each other by attractive 
Van der Waals forces. However, a principal method of analysis of problems of 
interfacial effects rests upon the assumption that a mean molecular field can describe the 
liquid. It is assumed possible to define an element of the liquid that is smal] compared to 
the range of intermolecular forces but large enough to contain a sufficient number of 
molecules. This approximation implies that on average the attractive force on any 
molecule in the liquid is the same in all directions giving to the liquid its fluid 
characteristics (Probstein, 1989). 

Mae Liquid-Gas Interface 

At a liquid-gas interface although the molecules are free to move in the liquid, 
their motion is far more restricted than in a gas where there 1s a little attraction between 
the molecules. Therefore the attraction between the liquid molecules will prevail and 
prevent the liquid molecules from escaping into the gas. As a result, the liquid molecules 
at the surface are attracted inward and normal to the liquid-gas interface, which us 
equivalent to the tendency of the surface to contract. 

If the interface is curved, a mechanical balance shows that there is a pressure 


difference across the interface. The pressure is higher on the concave side. Two surface 


tension forces balance the pressure increase in the interior. For an arbitrary curved 
interface whose principal radii of curvature are Ry and R> and have a surface tension o. 


A force balance normal to the surface will show that the pressure differential is 
1 1 
Ap= o{4 + z| (1) 


The previous equation is known as the Young-Laplace Equation. For a special case of a 
spherical bubble or drop in immiscible liquids, a well-known result is obtained, where r is 


the bubble or drop radius. 
Ap = —~ (2) 


When a drop is placed on a plane solid surface, it will be in contact not only with 
the surface but also with a gas such as air. The liquid may spread freely over the surface, 
or it may remain as a drop with a specific angle of contact with the solid surface. A force 
balance would involve both surface tension (c) and contact angle (@). If the contact angle 
is less than 90°, the liquid is said to wet the solid; if the contact angle is greater than 90°, 


the liquid is nonwetting (White, 1994). 
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Figure 1. Static Contact Angle Between Liquid and Gas Interface. 


ae Capillary Number (Ca) 

A dimensionless quantity which plays a significant role in determining the effect 
at the interface of the liquid-gas boundary is the capillary number (Ca). The capillary 
number measures the ratio of the viscous force to surface tension force and is defined by 


A small capillary number indicates that the viscous forces are not dominating the flow, 


ie viscous force HU (3) 


surface tension force oO 


the surface tension acts as a stabilizing factor in examining the effects of viscous 


fingering. 


B. VISCOUS FINGERING 

1. Theoretical Background 

An important instability phenomenon that has a widespread TiDDEenee in this 
experiment is viscous fingering. This is essentially an instability of bounded creeping 
flows, such as appear in a porous medium or a Hele-Shaw cell, when one fluid is 
displaced by a second fluid of lower viscosity. Differential pressures imposed at the open 
boundaries could be the driving mechanism for the flow. The fundamentals of creeping 
flows and Hele-Shaw flows are first briefly reviewed before applying them to the viscous 
fingering relevant to this study. 

Pa Creeping Flows 

Some approximate solutions of the Navier-Stokes equations for a limited case 
when the viscous forces are considerably greater than the inertia forces can be made. 


Since the inertia forces are proportional to the square of the velocity whereas the viscous 


forces are only proportional to its first power, it is easy to appreciate mathematically a 
flow for which viscous forces are dominant. Under these conditions the Reynolds 
number is very small and the inertia terms can be simply omitted from the equations of 
motion as a first order approximation. It can be seen that the incompressible form of the 


Navier-Stokes equations given by 


Dw 
PH, =~ eradp + BV Ow (4) 
with the inertia terms neglected become 
grad p= uv °w (5) 


The same boundary conditions as the full Navier-Stokes equations must be introduced to 
the system of equations, namely those expressing the absence of fluid slip at the walls. 

5 Hele-Shaw Flows 

A solution of the three-dimensional equations of creeping motion can be obtained 
for the case of flow between two parallel flat walls separated by a small distance 2h. Ifa 
cylindrical body of arbitrary cross-section is inserted between the two plates at nght 
angles so that it completely fills the space between them, the resulting streamlines is 
identical with that in potential flows about the same shape. It can be shown that the 
solution for creeping motion then possesses the same streamlines as the corresponding 


potential flow (Schlichting, 1979). 
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Figure 2. Hele-Shaw Cell. 


We select a system of coordinates with its origin in the center between the two 
plates, and make the x, y-plane parallel to the plates, the z-axis being perpendicular to 
them as in Figure 2. The body is assumed to be placed in a stream of velocity U, parallel 
to the x-axis. 

From the continuity equation, assuming two-dimensional, incompressible, fully 


developed flow 


In addition, the x-momentum equation simplifies to 


oP Oru 
ae . 


Integrating with respect to z and applying the boundary conditions, du/dz = 0 @ 


z=0andu=0 @z =+h/72 leads to a velocity profile 


eile a (8) 
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From the velocity profile the average velocity is obtained 
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4. Viscous Fingering 

The Saffman—Taylor instability at an interface could arise when two fluids of 
different viscosity are pushed by a pressure gradient through a Hele-Shaw cell. The 
instability is in many ways similar to the Rayleigh-Taylor instability: the primary 
difference is that the equilibrium state in the former is a dynamic one, in which the 


interface between the two fluids is moving rather than stationary (Faber, 1997). 





Nm Interface 


Figure 3. Side View of Radial Hele-Shaw Cell. 


Here we closely follow and adapt the development by Faber (1997). Suppose the 
cell is horizontal, consisting of two parallel radial disks, gravitational effects can be 
ignored. There is a pressure gradient, which drives the fluid in the radial (+r) direction 
with some uniform velocity (U,). At equilibrium, the interface between the two fluids is 


a concentric circle R = U,t. Where r < U,;t, the fluid viscosity is denoted by pw’; where r 


>U, t, the viscosity is given by p. Using the Navier-Stokes equations in cylindrical 


coordinates, the pressure gradient needed to maintain this motion of the two fluid regions 








are given by 
Op' I2y'U, op AG) (10) 
a ame or a 


Where U, is the average interface velocity and d is the disk separation distance. The 
pressures p’ and p are not necessarily equal at the interface, because the interface may be 


curved in the z direction. 


12uU 12n'U 
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Now suppose that the interface is perturbed in a way that at time t the nominally circular 


interface lies at r= R, where 


R os U i- aor (12) 





Interface 
Figure 4. Plan View of Radial Hele-Shaw Cell. 


There must be some corresponding perturbation in p’ and p, and it must have the 
same periodicity in the r direction. Any perturbation can be expressed as an infinite sum 
of wavelike functions. Since the perturbation cannot effect the pressure at large distances, 


these pressure terms have the form 








12uU 
p=- e - (r - U,t) ae (DF Ce 
(13) 
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where k = 27/i, and C’ and C are constants to be determined by boundary conditions at 
the interface. The boundary conditions, applicable in each case at r= R, and linearized by 


omission of the higher order terms yield 
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The pressure differential due to the interfacial surface tension is _- 
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The corresponding growth rate s, can be found to be 





ae, l aaa 
— ene aa - U, ku - “’) (18) 


It may be noted that the interface is stable for all values of k when py’ < up. When 
u <p’, 1.e. a viscous fluid is being displaced by a less viscous one, the interface can 
become unstable. The neutral stability curve corresponding to zero growth rate shows 


that it is marginally stable for the critical wavenumber k = k, where 


_ 12U,(u' ~ u) 
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(19) 


So for wavelengths oo > 1 > A, the interface is unstable. The perturbation that grows the 


fastest is when s, 1s a maximum, which yields the most dangerous wavelength as 





(20) 
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This corresponds to a wavelength (Amax) which is the fastest growing disturbance and 


hence the one most likely to develop into fingers 


he, (22) 


1] 


The largest value of the wavelength, which is consistent with the circular geometry of the 


flow cell of radius, r, is 


See = 2ar (23) 


Then for a given combination of interface velocity, gap separation, surface tension and 
fluid viscosities, the most dangerous wavelength can be found and the number of fingers 


that consequently appear can be calculated from 





i (24) 
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Figure 5. Growth Rate vs Wavelength (Radial velocity and Displacement Constant). 


Figure 5 is a typical plot of the growth rate as a function of wavelength for a 
given separation velocity. The plot clearly shows that the viscosity of the fluid affects 
growth rate and its magnitude increases with increasing viscosity. It can also be said that 
the maximum growth rate is not only larger, but the peak is also more pronounced, for 


higher viscosities. 


C. PREVIOUS STUDIES 

The paper by Paterson (1981) is the only earlier study to have dealt with radially 
driven flow in a circular Hele-Shaw cell. An approximate equation for the growth of 
viscous fingers was proposed. To model the initial growth of fingers, Darcy’s Law for 
governing the velocity of flow in a porous medium or Hele-Shaw was used. An equation 
for critical wavelength was obtained, and it was shown that if the circumference of the 
injected bubble was less than the critical wavelength the interface was stable and no 
viscous fingering was present. In the review by Saffman (1991) the selection mechanism 
and stability of fingers and bubbles in the Hele-Shaw Cells were discussed. The main 
focus was to better understand the uncertainties with respect to anomalously thin fingers 
and the nature of the finite amplitude of the instabilities. He suggested that the 
uncertainty in the boundary conditions at the moving interface may be the cause for the 
difference in the theoretical and experimental results. It was also noted that for unsteady 
fingering the problem is not yet resolved and many open questions remain especially 
since unsteady flows form singularities 1n finite time. 

In a recent experimental study by Maxworthy (1989), the study of the stability of 
the circular interface in a Hele-Shaw when air displaces a viscous 01] has been used to 
estimate the critical wavelength for large capillary numbers (Ca). A number of questions 
have arisen about determining the exact boundary conditions to be applied at the 
interface, which was slightly modified in this experimental study. The experimental 
results had a tendency to approach a wavelength that is a constant multiple of the gap 
width for values of Ca greater than unity. The reasons for this discrepancy are due to the 


non-zero finite nature of the Reynolds number, and three-dimensional effects. 
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Ill. OBJECTIVES OF THE PRESENT STUDY 


The objective of the present study is to obtain some fundamental insight into the 
onset of viscous fingering subjected to tensile forces. Unlike previous studies in this 
field, the main focus is the mechanics of separation of a thin liquid layer trapped between 
two parallel surfaces, which are being pulled apart at a constant velocity normal to the 
liquid layer. There is no shear or compressive forces acting on the sample. A 
mathematical equation is derived assuming a Hele-Shaw relationship for the onset of 
viscous fingering, known as the critical wavelength. A similar relationship will also be 
developed for the maximum growth rate of the viscous fingers, known as the maximum 
wavelength, which will then determine the number of fingers produced. The velocity of 
the interface and separation of the disks will be investigated. Experiments have been 
conducted and will be correlated to the mathematical equation to the onset and maximum 
growth of the viscous fingers. The work is experimental in nature and is aimed to gather 
both qualitative and quantitative data. Experiments have been planned and conducted in a 
way that some fundamental insights about the effects of critical wavelength interface 
velocity, separation distance and fluid viscosity. This study also aims to analyze the 


different competing regimes in the fluid deformation process. 
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IV. EXPERIMENT 


A. SELECTION OF THE MODEL SYSTEM 

Since this experimental study involves a liquid layer in tension (normal to the 
plane of the layer), the experimental apparatus was constructed to apply a purely tensile 
force in the vertical direction. The initial separation of the plexi-glass disks in which the 
liquid layer resides is about 100 pm. The stage that contained the plexi-glass disks was 
constructed to ensure that the mating surfaces of the disks were parallel to one another. 
The plexi-glass disks were also developed to ensure parallelism of their mating own 
surfaces. The vertical travel of the positioning actuator was aligned to prevent any shear 
force to develop during the experimental procedure. The alignment of the vertical 
actuator and plexi-glass disks were perfectly perpendicular to the direction of the ~ 
separation movement. Figure 6 shows the experimental apparatus, the supporting 
structure was constructed of high grade aluminum and mounted on a permanent rigid 


support platform to prevent any deviation from the normal. 





Figure 6. Experimental Set-Up, Apparatus. 
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The experimental set-up was designed to gather both numerical and visual data. 
One CCD camera was placed in the set-up to provide this required data. Since one of the 
objectives of the experiment was to determine onset of viscous fingering, the CCD 
camera was positioned to record the movements of the liquid layer from the bottom of the 
plexi-glass disks. This camera was mounted on the vertical moving stage (to move 
together with the disks) to prevent the focal length from changing during the operation of 
the experiment (Figure 7). The CCD camera was connected to monitors for viewing and 


recording. 





Figure 7. Movable Stage, Plexi-glass disks, Actuator, Microscope and Bottom View CCD Camera. 
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Figure 9. Experimental set-up showing monitor for bottom view. 


19 


The computer was used to run and store data for latter use and analysis (Figure 8). 
A software program written for this particular experiment controlled the vertical 
movement of the actuator. The software program enabled us to control the actuator by the 
computer and to maximize the accuracy of the experiment. From the control unit, 
displacement of the discs, tume and from the load cell, the force information was 
recorded. 

The experimental set-up was designed to give numerical data and visual data; the 
visual data could be viewed and recorded via the CCD camera (Figure 9). The numerical 
data was achieved through the use of a load cell, which was mounted on top of the upper 
plexi-glass disk (Figures 10 and 11). The load cell 1s connected to an amplifier and 
multi-meter. The load cell was used to measure the tensile load acting of the fluid; the 


data was then stored into the computer for further analysis and plotting. 





Figure 10. Close-up view of disks, microscope, load cell and lighting. 
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Figure 11. Close-up view from the top. 


B. EXPERIMENTAL SET-UP 

This section introduces each component of the experimental apparatus and 
includes a brief description of the component. 

1. Silicone Oil 

The “Series 200” range of silicone oils from Dow Corning was chosen to be the 
most appropriate for the purposes of this experimental study. They were commercially 
able to obtain, that provided a range of medium to high kinematic viscosities (50 — 
100,000 cSt) without any appreciable change in surface tension and density. They were 
easy and safe to work with, and they had advantageous chemical properties. Detailed 
information about these liquids is presented in Appendix-A. 

2. Positioning System 

In order for surface tension to dominate viscous forces, the Capillary number 
(Ca) must be small. For medium to high viscosity silicone oils, a small Capillary 
number (107 =i Oe) requires velocities to be small as well. Figure 12 presents velocities 
and Capillary numbers that are associated with different viscosity fluids. 

The positioning system control unit was chosen to provide these small velocities. 


The detailed information about the positioning system is presented in Appendix-B. 
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Capillary # vs. Velocity ( for diff. viscosities) 
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Figure 12. Capillary number vs. Velocity plot. 


3: Load Cell 

From previous work conducted by (Isik, 1999), the range of forces was already 
known. Concerns that were critical in determining the type of load cell were, sensitivity, 
minimum threshold level, resolution level, possible drifting in measurements with time, 
signal output, size, accuracy, and ability of being used in tension. The specifications and 


dimensions of the load cells are presented in Appendix-C. 
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4. Optic Systems 
Optical system of the experimental set-up is composed of a microscope that has a 
lens magnification of 50X and one CCD camera connected to a monitor and/or a VCR. 


The CCD camera was positioned to achieve a bottom view of the interfacial liquid layer. 


C. EXPERIMENTAL PROCEDURE 

Er ocrimente have been conducted for 50, 100, 200, 10,000, 60,000, and 100,000 
cSt silicone oils. For each silicone oil viscosity a large range of values were used for the 
initial liquid layer thickness and the separation velocity. A brief list of the variables and 
ranges are listed below. 


LIQUID : Silicone oil, constant surface tension & density 
LIQUID VISCOSITY (wo) : 50, 100, 200, 10,000, 60,000, 100,000 cSt 
LIQUID LAYER THICKNESS (1,) : 25, 50, 100 (125), 200 (250), 312.5, 500 um 
SEPARATION VELOCITY (V) : 1.2 um/s — 45.5 m/s 


Each experiment was conducted according to the following procedure and 


schematic diagram of the apparatus is provide in Figure 13: 


1) Experimental set-up was turned on to warm-up. 

2) Computer program was started and control was transferred to the 
keyboard. 

3) The plexi-glass disks were separated to introduce silicone oil. 

4) The silicone oil was spread on the lower disk by an injector. 
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5) The actuator was moved to position the plexi-glass disks to the desired 
initial liquid layer thickness. It was made sure that there was no air bubbles in the 
liquid, and that the liquid had spread completely on both disk surfaces. 

6) The velocity, step size and iterations were inputted into the computer 
program. The computer recorded and displayed step number, time, displacement 
and force in volts as output. Visual data could be monitored or recorded through 
the CCD camera. 

7) Each run was carried through until separation of the liquid layer occurred. 
8) These data have been analyzed by using plots of the force-displacement 


data force-Capillary number and the visual recorded data. 
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Figure 13. Schematic Diagram of the Experimental Set-up. 
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V. RESULTS AND DISCUSSION 


A. INTRODUCTION 

The experiments were carried out for a wide range of viscosities. A wide range of 
separation velocities and initial liquid ie thicknesses were also tested. Viscosities 
ranged from 50 cSt to 100,000 cSt, separation velocities from 1 to 47.5 tm/s, and initial 
liquid layer thicknesses from 25 to 500 um were used. Appendix D presents the 


combinations of these variables. 


B. FORCE-DISPLACEMENT PLOTS 

The following plots yield the resultant force as the liquid layer is being separated 
at a constant velocity. The force increases almost linearly to a peak value (Fmax) and 
occurs at a plate separation distance (dmax) in the liquid layer as explained earlier. The 
force then drops rapidly and almost instantly to a near zero (steady-state) value and 
remains there with no further appreciable change. Visually, the sudden decrease in force 
coincides in the onset of viscous fingering. The medium viscosity fluids (50 — 200 cSt) 
did not show viscous fingering to the same extent of the high viscosity fluids (10,000 — 
100,000), so the decline after the peak ae was more subtle. 

After reviewing the force-displacement plots, Fmax and dmax were dependent on 
the initial liquid layer thickness (1,), the separation velocity (V), the liquid viscosity (v), 
and of course surface tension (a special feature of these silicone oils is that the surface 
tensions is relatively constant throughout the viscosities). Each one of these variables 
were varied in order to achieve a relationship for Fmax and dmax. The following figures 
show the dependence of each of the variables for both the medium and high viscosity 


silicone ouls. 


27 


1) Fmax, dmax < U 

For a given initial liquid layer thickness (1,) and separation velocity (V), peak 
force (Fmax) and the location of the peak force (dmax), increase with increasing liquid 
viscosity (bv) as Shown in Figures 14 and 15. The magnitude of the force differs greatly 
between the high and medium viscosity silicone oils. This result may be interpreted as 


the extent of the onset of viscous fingering in the fluid yields greater forces. 


—e — Viscosity = 50 cSt 
—#-— Viscosity = 100 cSt 
—4— Viscosity = 200 cSt 


force (Ibf) 





displacement (mm) 
Figure 14. Force-Displacement plot that shows the effect of medium viscosity fluids on F,,, and daa, The 
Initial Liquid Layer Thickness (1,) is 50pm, and Separation Velocity (V) is 24um/s. 
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Viscosity = 10,000 cSt 


Viscosity = 60,000 cSt 


Viscosity = 100,000 cSt 


force (Ibf) 





displacement (mm) 
Figure 15. Force-Displacement plot that shows the effect of high viscosity fluids on F,,,, and d,. The 
Initial Liquid Layer Thickness (1,) is 312m, and Separation Velocity (V) is 12um/s. 


2) Fmax, dmax < V 
For a given liquid viscosity (vb) and initial liquid layer thickness (1,), peak force 
(Fmax) and the location of the peak force (dmax), increase with increasing separation 


velocity (Figures 16 and 17). 
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V = 22.7 um/s 
V=4.5 um/s 
V=2.2 um/s 


force (Ibf) 
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Figure 16. Force-Displacement plot that shows the effect of Separation Velocity on Fy, and d,,, The 
Initial Liquid Layer Thickness (1,) is 50m and Viscosity (v) is 50 cSt. 
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Figure 17. Force-Displacement plot that shows the effect of Separation Velocity on Fyax and dna The 
Initial Liquid Layer Thickness (1,) is 500um and Viscosity (v) is 100,000 cSt. 


3) Fmax, dmax & I/Ip 
For a given liquid viscosity (ob) and separation velocity (V), peak force (Fmax) and 


the location of the peak force (dmax), increase with decreasing initial liquid layer thickness 


(l,) (Figures 18 and 19). 
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Figure 18. Force-Displacement plot that shows the effect of Initial Liquid Layer Thickness on Fyax and dinax 
Separation Velocity (V) is 24um/s and Viscosity (v) is 50 cSt. 
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Figure 19. Force-Displacement plot that shows the effect of Initial Liquid Layer Thickness on F,,,, and dmax 
Separation Velocity (V) is 4.5um/s and Viscosity (v) is 100,000 cSt. 


4) A Physical Model 

A model of the force-displacement behavior in the separation mechanics of the 
layer can be developed based on a creeping flow analysis of the governing Navier-Stokes 
equations in axisymmetric cylindrical coordinates. The slow motion allows a quasi- 
steady analysis in which viscous effects are balanced by pressure gradients while inertia 
effects are negligible. The analysis follows the model problem suggested by Batchelor 
(1967) or Landau & Lifshitz (1987) and indicates that the separating motion of the two 
disks eas a negative suction pressure in the liquid layer gap thus resulting a resisting 
clamping effect opposite to the direction of separation. This resistive force 1s proportional 
to liquid viscosity and separation velocity, and inversely proportional to the cube of the 


liquid layer thickness at any instant, a feature that is qualitatively borne out by the current 


study. However a more thorough corroboration between model and experimental data 


will be considered in future work once more extensive data has been gathered. 


C. VISCOUS FINGERING 

It was proven that at the instant viscous fingering occurs, the interface in no 
longer stable. The interface motion is driven by the separation of the two disks. Surface 
ies serves as a Stabilizing entity between the two immiscible fluids. From the 
discussion in the early chapter on the Saffman-Taylor instability, a critical wavelength 
was determined at the onset of viscous fingering. At the onset of viscous fingering a 
number of disturbances develop along the interface called corrugations. At this point 
the surface tension can no longer stabilize the interface and it is no longer undistorted and 
circular. The interface is now unstable and continues to grow unstable. If the instability 
proceeds the interface is no longer balancing the difference in viscosities and pressure 
differentials to surface tension, the number of corrugations has reached its maximum 
value and coincides with the maximum resultant force. “From the Saffman-Taylor 
instability analysis of the experimental radial Hele-Shaw the maximum wavelength and 
number of corrugations were derived. The maximum wavelength is the largest 
wavelength that gives rise to the viscous fingers that grow the fastest. In order to prevent 
viscous fingering from occurring and leading to an unstable interface the critical 
wavelength critical wavelength becomes an important measuring device of instability. 

A comparison of various separation velocities and initial liquid layer thickness are 
considered below. The maximum wavelength and number of corrugations for each run 
was calculated and compared to the picture at the same instant. Table | contains the 


calculated quantities. 
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Kinematic 10000 10000 10000 10000 
Viscosity (cSt 
0.021 0.021 0.021 0.021 
/m 


0.0254 0.0254 0.0254 0.0254 


Initial Liquid Layer 0.0002 0.0002 0.0003125 0.0003125 
Thickness (m 






0.029 0.051 0.0328 0.051 


m/s 
m/s 
Wavelength (m 
Pa oe 
Corrugations 
Table 1. Experiments conducted for 10,000 cSt Silicone Oil, Calculating Maximum Wavelength and 
Number of Corrugations. 


As initial liquid layer thickness (I,) decreases, viscous fingering becomes more 
pronounced, leading to an increase in number of corrugations and a decrease in maximum 
wavelength. Similarly, separation velocity (V) increases, viscous fingering becomes 
more pronounced resulting in an increase in the number of corrugations and a decrease in 
maximum wavelength. Also, separation velocity (V) and the liquid viscosity (o) also 
affect the degree of viscous fingering. Higher the separation velocity and liquid 
viscosity, greater is the extent of viscous fingering. 

The following pictures were taken at the maximum growth rate, although the 
actual number of corrugations can be difficult to identify the same trend as shown above 


for the maximum wavelength and number of corrugations can be seen. 
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Figure 20. 10,000 cSt Silicone Oil, Initial Liquid Layer Thickness of 200um and Separation Velocity of 
12um/s (Scale of 1div=1mm), N80. 


Figure 21. 10,000 cSt Silicone Oil, Initial Liquid Layer Thickness of 200um and Separation Velocity of 
24um/s (Scale of 1div=1mm), Nx84. 
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Figure 22. 10,000 cSt Silicone Oil, Initial Liquid Layer Thickness of 312.5ym and Separation Velocity of 
12um/s (Scale of 1div=lmm), N=47. 


Figure 23. 10,000 cSt Silicone Oil, Initial Liquid Layer Thickness of 312.5um and Separation Velocity of 
24um/s (Scale of 1div=l1mm), N52. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


Well-defined trends in the force-displacement behavior of thin liquid layers were 
identified from the analysis of the recorded quantitative and qualitative data in this 
experimental study. A wide range of viscosities were examined and all force- 
displacement curves have lead to the observation of the presence of a peak force (Fmax) 
whose location and value depend on initial liquid layer thickness, separation velocity, 
liquid viscosity, and surface tension. 

The presence of the peak force is more pronounced and its occurrence coincides 
with the onset of the viscous fingering mechanism. The approach to the peak force (Fmax) 
on the foreeuaieelceniea plot, is accompanied by an instantaneous and sudden decrease 
in force which is followed by a near zero (steady-state) value with no further appreciable 
change. The linear stability analysis allows a fairly accurate prediction of the onset of 
the viscous fingering behavior and correlates well with the experimental data. 

Further research on this subject is highly desirable to widen the parametric study 
to other viscosity ranges and surface tension. A more accurate characterization of the 
interface velocity is crucial to further the understanding of the onset of viscous fingering 
and its dependence on the Capillary number. A review and refinement of the analytical 
model may also be conducted and balanced with the necessary boundary conditions, 
which could lead to corroboration with the different observed regimes. 

In addition, a more advanced and integrated testing apparatus with better 
measurement abilities of the interface velocity, layer thickness, actuator motion and 
interface flatness and parallelism is required. It would also be desirable to increase the 


data acquisition rate to at least the milli-second level in order obtain more accurate force- 


of 


displacement information. The use of an optical grade glass, instead of plexi-glass, for 
the disk would be advantageous to ensure smooth, flat and parallel, and hard (scratch- 


resistant) surfaces. 
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APPENDIX —A 


SILICONE OIL 


Medium Viscosity Silicone Fluids: 

200° Fluid, 50 cSt 

200® Fluid, 100 cSt 

200° Fluid, 200 cSt 

High Viscosity Silicone Fluids: 

200° Fluid, 10,000 cSt 

200° Fluid, 60,000 cSt 

200° Fluid, 100,000 cSt 

Description: 

200° Fluids from Dow Corning, 50 — 1,000 centistokes (cSt) are medium, and 
10,000 — 100,000 centistokes (cSt), are high viscosity polydimethylsiloxane polymers 
manufactured to yield linear polymers with average kinematic viscosities ranging from 
50 — 1000 cSt and 10,000 — 100,000 cSt. 

Composition: 

Linear polydimethylsiloxane polymers characteristically have the following 
typical chemical composition: 

(CH;) SiO [SiO(CH;) 3] ,Si(CH;) 3 

Commercial bulk polymerized dimethl silicone fluids, such as 200® Fluids typically 
contain trace amounts of process impurities. 


Benefits: 


200® Fluids, 50 — 100,000 cSt, have the following product characteristics; 


> Clear > Nongreasy > Nonocclusive 


> Non-stinging on skin 
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200® Fluids, 50 ~ 100,000 cSt, when compared with other materials that may be 


substituted in a given application, may offer one or more of these comparative 


characteristics: 

> High compressibility > High damping action 

> High oxidation resistance > High shear-ability without breakdown 
> High temperature serviceability = High compatibility 

> High spreadibility > High water repellency 

> Low fire hazard > Low odor 

> Low reactivity > Low surface energy 


> [ow temperature serviceability = Low vapor pressure 
> Good heat stability > Soft feel and lubricity on skin 


> Good leveling and easy rubout 
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APPENDIX -— B 


POSITIONING SYSTEM 

Model 855 Programmable Controller System (Newport Corporation) 

The Model 855 Programmable Controller is a microprocessor-based system that 
allows simultaneous direct or programmable control of up to four Newport linear 
actuators or rotary stages. It is simple but powerful vocabulary of mnemonic commands 
allows straightforward programming and control via the handheld 855K Keypad/Display. 

The 855C Controller: 

The 855C Programmable Controller is the nucleus of a system that automatically 
controls up to four Newport precision positioners and stages. Its large, easy-to-learn 
instruction set and standard RS-232C and IEEE-488 interface ports allow it to work 
closely with external computers and other data devices. Its programmability provides 
stand-alone automatic control of actuator motion, yet no knowledge of programming 
techniques is required. It also supports the optional 855K handheld Keypad/Display for 
convenient data entry, control and program editing. 

Specifications: 

Actuator Control: Provides +15 VDC, 1.5 Amp power and +10V velocity control 
signals for up to 4 Newport linear positioning devices. Receives 


and decodes dual output, 90° phase, +12.5, -0.5 V encoder pulses 
and limit signals. 


CEU: 8 bit 6809 high-performance micro-processor 

RAM: 4 Kbytes 

PROM: 24 Kbytes 

EPROM: 2 Kbytes 

V/O: RS-232C Serial Port/IEEE-488 Parallel Port 

Power: 110-220 VAC, 50/60 Hz and 850 mA with all actuators fully 
loaded. 

Dimensions: 5.25” high X 17” wide X 12” deep 
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The 855K Keyboard/Display: 

The 855K can be used to completely control and program the 855C. It is handy 
for manual control of actuator position and for the reviewing and editing 855C programs. 
Its backlit alphanumeric liquid-crystal display shows position information for the four 
drivers, programming instructions for the 855C, and messages from the 855C. The 
855K’s keypad has 35 durable, tactile-response membrane keys for command and data 
entry. The 855K connects to the 855C via a coiled cable. 

The Linear Actuators: 

The 855C is compatible with Newport’s 850 Series Linear Actuators. Each of 
these devices use DC motor drives with integral optical encoders for smooth operation 
and high resolution. The 855C fully supports their resolution and range. Only one 


actuator is used for this experimental study. 
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APPENDIX —C 
LOAD CELL 
Since the viscosities of the fluids differed greatly from 50 cSt to 100,000 cSt, 
there was a need to use two different load cells for the experiments. For higher 
viscosities an MLP-100 (100 Ib maximum load) and for the lower viscosities the MDP- 
2.5 (2.5 lb maximum load) were used. The specifications are listed below for each load 


cell. 


MLP-100 Specifications: 


MINT LOW PROFILE LOAD CELL 
UNIVERSAL | TENSION OR COMPRESSION 


The model MLP Series load cells were designed with economy as first priority. 





The MLP Series are anodized aluminum with a unique low profile design, which 


provides excellent stability for in line application for tension and/or compression. 





10" & CONDUCTOR COLOR COCED, 


SHIELVED CABLE SUPPLIED 
THREAD TYPICAL 
BOTri ENDS 
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SPECIFICATIONS 


Rated Output (R.O.): 2 mV/V nominal E. 
Nonlinearity: 0.05% of R.O. | ae 
Hysteresis: 0.05% of R.0. s 
Nonrepeatability: 0.05% of R.O. ig 
~ Zero Balance: 1.0% of R.O. — 
Compensated Temp. Range: 60° to 160°F w 
Safe Temp. Range: -85" to 200°F = 
Temp. Effect on Output: 0.005% of Load/F ig 
Temp. Effecton Zero: 0.005% of R.O/F aa 
Terminal Resistance: 350 ohms nominai E 
Excitation Voltage: 10 VDC 

Safe Overload: 150% of R.O. 

Weight: 107. all ranges 





iy pap ae Pl oR Ar ee, Po Mesa) eee on “~ 


as) # J 
1 hiah tsi I OPT CR COMA OR ER EMOTE ERG BOER SE ae HO se yioesy 3 
ESS py SEAS SENN AYR ae RA ie = So ad oD, OY Ee an, es od TESS is 
1 EN PAR ee Penny Oe EIA aici piss af Bee ie Tar alo las AIS ees yo pet ae ae 


DIMENSIONS (INCHES) 


maior eo ee =? 


|MLP-150 | 150 2 : 
i200 “ies so st ome seat 
300: z 
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MDB-2.5 Specifications: 


ULTRA PRECISION MINI LOAD CELL 
UNIVERSAL } TENSION OR COMPRESSION 





The MDB Series was designed to help fill the growing need for a greater selection 
of high accuracy load cells for use in space limited applications. The anodized aluminum 


MDB's are compliant tension and compression, therefore, a good choice for in line 


through zero applications, as well as single direction tension or compression. 


MICROTECH CONNECTOR 
WH 10° 4 CONDUCTOR, 
COLOR CODED, SHIELDED 


CAGLE SUPPLIED ae 






pes 





\ v4-28 UNF 
220 PEEP BOTH ENDS 


See desc 


"eerpernraty EE GME te nee 


~ Rated Output (R.0): 
| _ Nonlinearity: 


| Compensated Temp. Range: 


‘Sate Temp. Range: 
Temp. Effect on Output: 


“Temp. Effect on Zero: 


~ Terminal Resistance: 


Excitation VoRage: 
Safe Overload: 


Weight: 


“2mVN nominal 


0.05% of F.0. 


: 0.05% of R.0. 
> 0.05% of R.O.. 


1.0% of R.0. | 

60° to 160°F 

65° to 200°F 
0.005% of Load/F 
0.005% of ROSF 
350 ohms nominal . 
10 VDC 


150% of R.O. 


1 02. all ranges 
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APPENDIX — D 


EXPERIMENTAL DATA 


Initial Liquid Layer Velocity 
(uum/s) 


Thickness (I,) 
25 um 


Initial Liquid Layer Velocity 


Thickness (1,) (um/s) 
SOum 


Initial Liquid Layer Velocity 


50,100,200" | "a (um/s) 
cSt a 


Silicone Oil 


Initial Liquid Layer | Velocity 


Thickness (Io) (m/s) 
200um 


Initial Liquid Layer Velocity 


Thickness (],) (jum/s) 
312um 


Initial Liquid Layer Velocity 


Thickness (1) (m/s) 
500um 





Table 1. 50, 100, 200 cSt silicone oils and their varying parameters. 


47 


10,000, 
60,000, 
100,000 cSt 
Silicone Oil 


Initial Liquid Layer 
Thickness (ly) 
125m 


Initial Liquid Layer 
Thickness (ly) 
250um 


Initial Liquid Layer 
Thickness (1,) 
312um 


Initial Liquid Layer 
Thickness (Ip) 
500um 
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Velocity 
(u/s) 


Velocity 
(um/s) 


Velocity 
(um/s) 


Velocity 
(um/s) 





LIST OF REFERENCES 


Batchelor, G. K., An Introduction to Fluid Dynamics, Cambridge University Press, 1967. 
Landau, L. D. and Lifshitz, E. M., Fluid Mechanics, Pergamon Press, 1987. 
Faber, T. E., Fluid Dynamics for Physicists, Cambridge University Press, 1997. 


Maxworthy, T., Experimental Study of Interface Instability in a Hele-Shaw Cell, The 
American Physical Society (1989), vol. 39, 11. 


Paterson, Lincoln, Radial Fingering in a Hele-Shaw Cell, Journal of Fluid Mechanics 
(1981), vol. 113, 513-529. 


Probstein, Ronald F., Physicochemical Hydrodynamics, Butterworths, 1989. 


Saffman, P. G., Selection Mechanisms and Stability of Fingers and Bubbles in Hele-Shaw 
Cells, IMA Journal of Applied Mathematics (1991), 46, 137-145. 


Saffman, P. G. and Taylor, G. I., The Penetration of a Fluid into a Porous Medium or 
Hele-Shaw Cell Containing a More Viscous Liquid, Proc. R. Soc. Lond., 1958. 


Schlichting, Herman, Boundary Layer Theory, McGraw-Hill, Inc., 1979. 


Isik, Sefa, A Preliminary Experimental Study of the Behavior of Liquids under Tension, 
Master’s Thesis, Naval Postgraduate School, Monterey, California, June 1999. 


Sherman, Frederick S., Viscous Flow, McGraw-Hill, Inc., 1990. 


Taylor, G. I., Cavity Flows of Viscous Liquids in Narrow Spaces, Proc. 2° Symp., Naval 
Hydrodynamics, 1958. 


Taylor, G. I. and Saffman P. G., A Note on the Motion of Bubbles in a Hele-Shaw Cell 
and Porous Medium, Q. J. Mech. Appl. Math., 1958. 


White, F. M., Fluid Mechanics, McGraw-Hill, Inc.,3" Ed., 1994. 


49 


THIS PAGE INTENTIONALLY LEFT BLANK 


50 


INITIAL DISTRIBUTION LIST 


Defense Technical Information Center...................cccccc cece cess 


8725 John J. Kingman Rd., STE 0944 
Ft. Belvoir, VA 22060-6218 


Dudley Knox Eibrabye=...... 20232... .5) eee ee eee 


Naval Postgraduate School 
411 Dyer Rd. 
Monterey, CA 93943-5101 


(hawman, Code Miwa)... nants hee tcc. ee eee eee 


Department of Mechanical Engineering 
Naval Postgraduate School 
Monterey, CA 93943-5000 


Professor Ashok Gopinath, Code ME/GK....................0000 0 eee 


Department of Mechanical Engineering 
Naval Postgraduate School 
Monterey, CA 93943-5000 


Naval/Mechanical Engineering Curricular Office, Code 34........ 


Naval Postgraduate School 
Monterey, CA 93943-5000 


Commandme Otncer (Code 35) 5. eee ee eee 


Naval School, Civil Engineer Corps Officers 
Naval Construction Battalion Center 
Port Hueneme, CA 93043 


Chraste phen ll skemKOpen2. oi... cil Sscabeis ce eee ee eee eae 


6104 Summit Pointe Road 
Harrisburg, PA 17111 


ol 


























a = “ae 











ee 





















= . . 
®-5 ve wi le 
Shane Cn eae - 4 
+. + » « ae . 
- sm ef - « 2 ee " arises : 
te Ce ae See we oe F Cyawiete 
an et . - at ©. « =e oe te Me OLP ity * 
Sr Pee bees oe wares Ce ee emate ety ¥ Cae re 
' "= A fo- vas . a eee ar | 
ar Meee 2 Nee) Pam Si heaniee 
rm wwe oM- +. wen ow 
aie a=, fur oe, ye > 
> re “oe qu 





















1 ’ Brel es 
- nS . 
“> ~% oe . . ’ e F < 
. Lich NINA | « Be 6 . . . ee . Fs 
SD eee Sa, va cua Sa hae “ ais Teh eg Oy a a ae 7 é or it 3 
fs ee Ye one -m , Ss sa . See ha ee aS aS Shek P 
* *haeene x . Siena Us 7” Car a 4 a t is a a oe e we rads . i 
< Fam eye 7. - +. Spies oa <The =: “re - Pr <i =e ek . Saas x Ks e aa : 
>. "as sober ee UP Sy See S <= wee? tien & Cat acy > a ay ¢ = Sas a ashe Paid 7 geet) 
“05 -p meme - es 8 = .. embe Sige wie aca, ~ ne eee Gg nt a arte ae nt. aa 2 Pas ene 
1M « Pee ene ee. ee ae oa Pred . SENN e. Wiaislors . . > Cart ae are 
2 we w+ << wWigns = = au clear Ee tra ae “8 a he ae ws =F 4 
pS gsm ey ors site im Las weEe ee . - a s r) 
~~ es Nene C, ee a er ees Ae « 
es + * Nae P Ree PRL®.) S58 gre wapire ie rs - 
P vteen we ° - mt vie 
«> 6 - 
































"=. 
Sea SO OMe e. na St wee 
eS ater ek 
ve ~ 
~ < ‘ - 
a PS ‘ o a Ri taiay oa 
ie Bee e<~ . - . ” 
aye 7. as . "ss - = ‘ e . ean as 
-. tee BONO 5) Ce ‘ < - 
: Rear : . . 7. ° é . " - ; : se ° * 
s ‘ aa “« - « ¢> 2 ao . S +e ee ee . 6P- siwg 
. , . Cae A Je Ok hele ai ° ° , « -¢ 7 = * Oe Moi” it 
eS a he ad a. o * . ° > * 3 nd : 
A m4 cS . ve a ae ee O80 ee renin 
. é . = a +. tee . + hae a . ’ ods ade eee . * ge - . aed o? ‘ 
SRP wneen oye wee = Sane iad ne ‘ - +» Cn a ae. 4 er ey eee es ‘ - * ‘ oo - cara 
. es eg Wt Jet Sy a ~ ™ . * “ses “ee Soar. a . wee Cae ees Neue oe ‘a . 
Pati tng on nn eee Py, e wgen tom oe Sines) aeaae 80a ED Sola ig . Fee os cigs Sica oh asb es faen ba gels P é if 
«© 8 ey ieee he Le . 7 wee ore wee 8 Memes agg tla), Se, et tas . . . oe he ous 
Toms ge. -« Wis Cr Fig) is -* Pe Se es “e- eh Me te, Flag ca, Pace he * « 
ren 7 TP ree 26 | Ne we aw z Se xs Fe age 5 + ¢ > 
ee aon, w OT ay my weenie, ara ~ *ef,e e roe sp 
-tme eee tte > a t ® = m+ 
= 2 Bs@rune em mee ae, ; ce iat ts 
Fess mn ps ~ do cies. ek eres 2 ate wy 
Ts meh mem, “A Pas “ 





~s 

at 
a) 

07 ge y 


vs 


















by ror 
. . » «< 
ae ¢. 4 “be wera - « 
« we . *. xa 2 Shia 
Ls . =p S * " aye . 
wee erm zee 4-, 
vay et iia en at See eee eee 9 Bee tia 
ae gi ee ay FES ON cig fo . Wey isiteg hore 
-« Faeneaing, - . eee = ° Bee 
BPO hema ¢ 
a a é 
C988 pin mn 8 abies 
Prin te ah Bwee, ayy ase ed tet wee 
© Ort Me . tates ~~ 5 se + 
= By bP oy re) Ha nanan tt» 
- — 8a wy ig wT Pepi « © a 
—— BP mae ne 
OO —— 
—EE 





| 


KNOX LIBRARY 
| 
\ 
0041064 


1 
| 









Hi | 
| 
HL 


| 





1] 


= =F oy. 
Pats Hie phe w+ 
vt we 








role ie 


| 


THIN 
| 
3 2768 


| 
| 
| 





om, 


tee me meng 





@ vertey ok 


ts, 


| 
| 
| 





| 





Foyer ae 










~~ 7e™ ant Pee 
we : 
pT A er sean 
~ . 888 
‘ 
SS 
————— 


Serial Sty “sak 
Meee etane 2 


| 





| 


Mi 


Se Per anee 


WHT I| 
Hi | 





| 





| 
| 






| 





ures 


| 
] 


Oren we 
Pr 20s seathe « 














sFeere, 





~ seem ae, 


Shee asan iy 


soemh wm er, 
ht ene 
Reza 





Be PE 
* pinot 





= SPatet s 
whe 





Fate 
* Bae wiasc * 












sete ee 
tear 
ePare eye op 


ae 











often mie => 


te? 


. 
ed TLTTY 
TAeee 





‘ee 8 
-* . © ae oo © fae oy . ‘ 
She ts yee 
, 


Pere 





eT 90 oi 











“8 AP ep eure, 

























ae 
“* sPatecas ¢ 
eboney 
+ mtg op - ate 
sits My une eyys Ce ee ees Sd ed +4 vee ae “oo 
Tene Pane Te even a> 4 oer y VePage o>. PF hr tye, 
PA Puree. 7 eee FoF esata. te pee pte 
Cr es BS. rer edie ure. ag 
een FY Oe, aoe, Bie S Sores satan 
CORMIER eye TEN Bh oe ad 
"Fee ede d 28 - 





AON Rete we , 


MUS Orde peony 
eh See 


SP NE em way 


“SE Sater en 











bah 4d Leelee 


















SOP ee em ae 
noha eee 

SFr aneeny 

es isi" Sar 3 = 

<n Perper ratte Meni rn gay, 

sern Sree OR Beye asy ye 

ay 

WARY Won 





Pe AM ans aye 

















are gr a, 
San Pia eem "he Per gy 
° tpt seorge 
oe, 
a ‘ 7am ft eo same rase ¢ 
TENG lem “erst ueD pee toe we CS oe 
t= of ° oS et my TeMsentezedueun og -ot 
Pemee om Sot eRe ote os “eee melas Feo Org, a Mea, 
ore Paty, Taree tEte? eee ca Wee rens Ace, 
RBar. Patel hd a PT 
heddee SL es 





= ept0aza.e.e & 





ony 
dae 








Te ae the ey 
bind ide a Se 





PePeSem slay 






oa 









































or 
Se Om tesa, © oan 
fwees, neee, weugee Piysengh 
Ae he eee et 
del bar LOL tee Dem 
EER PAR sie gg vey oe we 
bata ied thar omy, % 
card * 6% ose ue hime 
Be kerae  C a ¢ # OAstrames . ay 
bart etal. . 
sPease a ‘ 
75 Gmete 
= bel hd ' 
a a See D 
3h" SESE RS eee muse eae oie oe 
Heme wae ada te afarer oa. 
FOS WORM OSA Mb eee TMF ogee 
Sed eT ry bie eid ee eT ers 18.7 as o 
SP beer SLR EM ony . 
Celome ms ser Bet er. 
cat” Pa eng ome 








s 4 









a ree 





+74 bo eumag 
2 Nh pagers 
MEP hae 



















































































































































































































































































wet owe 
ue + 8 gt? Peete 
ones *AtNere ote 
ERAS SA cay BPR teernens a, o <8: 
Yorn ee 2408 © Ae, a a 5 a4 a 
‘ ar a 
rae a | ae Saket eten ¢. Boawe 
= ®e ites 5 Sass ‘ ° he 
Mpa eicenek cae ieee Sia ala Mey, 
Sees e aN bt ahtdg Rs deeeie, 
bid SA oh Pe “ 
bo ~ Piao Ps a 
le? eer a belgie adda Ms Tn wtarye 
“a a8 ‘ - A MSeke 
cht een porte eg : 1 abe ne wmanen soz pe “Se 
Se A@ pa saegre, J ed Bed : we ase, Teh rre Seeyae 
Bay OTNIOF PLR Sean Oem tre. Ree 7 - 
AML Pua, 2 Pwen Suan any ae 
siete hte oe ed Meee, * eab wap FMS e tense Wesevey ay 888 deers woe 
POPU! + wrerery 250.) bihaiineheea hia ete! ores w oye ate eg - eT ee Seow age Pie . ene en ate ater a gs 
windin ake et 2 Merce Z Fic Si pai . *% Se ae ° tae. ere ahs abe ogepepe sie8 2 
ete © 2 tBomayri rye Gree > Pe ware OPS ts 0 os seteebe. Scena, Sheaves Bene Se LATE IE i SNe ares 5 tange ars Sehinewae spose cnet Shee dents hal? tan gee we 
ese ¥9 SOF enn, ‘ ~¢9 ss t oo eae =e . an SS Fld “ cad Cope ef net Wie 
Seen ts © mitts . es Sade ce lene 8 er a ALM canes : fis rary Gas me aki ue foo Oe Seely ave SFenes , Syrar anes wr ohare 0 a NS mw NE bee tay a 
» et SN nee annsas im ote eke ek ae pees desert wa Oe mde ne . : ane a Wis gts” AM Aas m ay Baers 
YW Mh er.eee > Peie ff Sen we 8 Fe yg . s tea okt kee | Med hap = hae hee ted Soe © wen, Qomae 
- - SHA bat Be ada, < éo0 Yeneers an Nien 
tet bs ee - Satene as Vree.y ”, “ . 5 okt tee, 
Attia we atentis ve, ah er te tL MPS tavengens Se phere ot he. OR et wee, as at eee ase Sates oan 
Paty Te aes. 287 oe ned urate, law er oo te os at Si mapweneare Ae eee tal aa oSz uy es PALL in cy ee ee P Sa aK eye, “og Sa, 
ees ser arte Ate thane: + oe Pye. 18S abey ary = 08" vage Scr Bigs vie Bai eA sours Ss ~2 ah pe t+ Ram on 
WCC EE L® Bike mS. toe wae Gg tong L eM eye fa’Sj ane- 
Ss WS NewS aRL aloes ag SFM, Nd ae Fade we ote ible ooh it Par Tp afene “a Paryy wt sizeg ri ge, 
Fa Bee age oF! a ibe Lett eT ere hte Let) awe 4 Few ecnenagy We eating OT ry “eae n a ns, © Gene oteaie. rane De 
7h omaure, r temegnen oe * 4 2 bo tacmane g See ee ® 
: re 4 STewere 5, ¥ a Pn xe . ‘ te bey ane P%e% ean age 
Ft ve ge he shh abt JO, al a. Bee. 
ie * «fe PO ees oe 
tPer Pera ian y te 
Spear we “obe* ewgnes oweage aot ae 
Pe» e I a* 2 DB cups Pete 
tir nhs, OY Bor dagen, yg, Gites BL Deets spar CeEPe elanye payne 
Pet ons, FOP estan Ln EPataseni a, ae IE Miers =" det : UO SEN ele ares ner, helo 
a! . Pe ~ . oeekan, 5 Fs eee aL “8 Figs « ts ofla= - : 
wees hy H 8 ag inet ie hadedel Labatt é ECan oe : % A Ras i Tr Sees wink ae we ae we Pe whose eb aD 
PEPE E., bbl al eee Oe Pak ed ge, NR tess. é? 6 pen tAece a wie Oy 
- 4 “+ ae se deny 
PE - Nay Uae ae, TE AP. Eee Dae ‘ as on bate 70 e Nesp en “Zarate ye amy 7 Pay Bi ane ea 
of ee be Sliabtaatal a te et 5 TP eten, Sete tae 3 bese es 4 Asauen nt al SE, PCT « fits? wave 
aad ales hace ane ee ei sae Sarasa yee, wt ay ee, $i. eh bs on, ae ees OO Seema a oe i "s ree ors . ae ; 
, » we i: cee + : é " => et = tee “te 3, ae teole ug o 
SOS RATAN, Nae erL ty, teed eer Pte irre ered iat MZ chee ayer eee aeaparee “u Bisne we cubase ke aM a Kee 
Orsay Amobbchte BOF om as. gy Loy 2 Dare, “Oe esa od eee LIT aren * ares Fate 
Wd teh td | Mere Ostas themes SIs 8 aN Bag, ore) oe ree os 
we "Tat asetay Motes 4s, a" Fem dane ts. At tare st 8 wy re . 
Var usame, Leet * SFE dees Wire gs Saou eta: oh) 
ee, 2 Pros tbat, bie oe ee 
Pete ele ope hb he ae Pee oS eB te 
Le Me am sae WERT Ber rie fue ORO ah gaa Fong - Wee eteey 
Fea Hee eer eres ya ny 
Le Ltd 










eu 2 


© Opes 
Mee wh oe taay 
. 


Sates Sy ot 
‘o> #06 


Maes aun, 
Se ee 






a9e ¢ "mapa ® 
Oe 


bl Ate Peet LLY 








aury 
bia ee ES 
wey 





@ 4.5 
Sehess saose, 
ess 


“tk 6 








“AP WS gin 
mie be te ee 





- SWF qq ne, ke 
> a tat be Ld * . 
: °. ars 
We os rae 


Wee woe ara 
seh 
o@ senarnpe, og 








P are- 
ott 
seeeh oe 











7. 
“e"tean 






S9+ Pye 
~ BP ee any 
owes 
*Georsye 


adie” a ots By 
WO = Es vote a, 


Se Pevaetse B %,, Bat 
Se 


Wer 8 





ates 






POO aD me aes on dunes 
thes De er eeree ys att | 
ee abe eet Oe} 

PSAREP Rerahy 

HP hemm O are ON wns tei hiaye 
tay 88 dee eny Taw gmine; 
PAON EE re Bary 
Pome sa: an . 











"thew wr eo bane, 
* by ary 
" ? bes 

NGhagy 
vee 





) 

Press « 

* Place ~ 
i 


4 

at 

sed nte as 

ewe 
“tte keer oe, . 

IPP may 

«4 


eiipptaehitaa SOL IPT Ct 





SP ase> pe ay 
Pee 





eee tues 
ten ea ata? 
oh stat bas 










odie 
Pay Vu on. 


*e8 ge 
5f4 0 see, 
one 

















es OP fs +t oon y 
~\ SUN aay ca Talag-~ 


eal lak tend 


















gehts Stam as 

fart te oe, 
OP Dawe age 

*tCbis wy 

SObds, eRe’ 

mae OL 

PtOae, 








: 
© WM sy estae, 
WOOT EP Urteh ety ow uae 

tatee, 


afte hei e 
“~ 

















OF De Ote74, 50 6 
ela 











“Me ame 9 
Hey 





‘ thea 
*t &a 
PonpaPatien. oy 


e aie 
VORP He as mes ae 





Mak Sany 


> Wwhet@rtn® oy 











Pit ablebhddatncnl of Laeted 
+ Seta 


Pt. 








